Universal threshold enhancement 
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By assuming certain analytic properties of the propagator, it is shown that universal features of 
the spectral function including threshold enhancement arise if a pole describing a particle at high 
temperature approaches in the complex energy plane the threshold position of its two-body decay 
with the variation of T. The case is considered, when one can disregard any other decay processes. 
The quality of the proposed description is demonstrated by comparing it with the detailed large N 
solution of the linear a model around the pole-threshold coincidence. 
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The near-threshold two-body decay of a resonance is 
a well-studied problem in scattering theory pj. In a 
medium with tunable parameters (e.g. temperature, 
density etc.) new features may arise @. It has been 
found 0, that a hallmark of partial chiral symmetry 
restoration in a medium is the enhancement of the scalar- 
isoscalar spectral function near the two-pion threshold. 
In order to elucidate further the origins and the nature 
of this phenomenon, in this note we formulate a general 
approach to it. On the basis of this investigation some 
features of threshold enhancement appear to be univer- 
sal, showing up whenever the mass of a stable particle can 
be driven by tuning some parameters of its environment 
to pass the threshold of a two-body decay into stable fi- 
nal particles. Following our general discussion we shall 
show evidence for this situation in the linear a model by 
studying the propagator of the a particle, neglecting the 
finite in- medium lifetime of the pions. 

Let us start our general discussion by assuming the ex- 
istence of a branch point in the complex ^Jp 2 ~ rest frame 
energy plane on the real energy axis for the propagator 
G CT (p 2 ,|p|), {p 2 — Pq — p 2 ), which corresponds to the 
a — ► a\ + a 2 hypothetical decay. (In this part of our 
discussion a refers to an arbitrary particle. The extra 
dependence on the 3-momentum is a medium effect.) We 
can restrict the discussion to the halfplane Ke^/jP > 0, 
due to the symmetric behavior of the propagator. At 
the branch point a square root singularity, proportional 
to 6 = \/T— p 2 /M 2 , shows up in the self-energy con- 
tribution to the propagator, where M denotes the tem- 
perature dependent position of the branch point. The 
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analytic (Riemann-sheet) structure of the -y/p^-plane is 
specified by requiring that G a is analytically continued 
through the part of the real axis where p 2 > M 2 , i.e. 
where the complex threshold factor S is purely imaginary. 
Furthermore the plane we are working on is bounded 
by a cut along the section p 2 < M 2 of the real axis. 
Along the cut 8 is real positive (negative) on the upper 
(lower) halfplane. Conventionally, one associates the up- 
per halfplane with the first, the lower one with the second 
Riemann-sheet. We shall call the upper edge of the real 
axis reached from the first Ricmann sheet as physical, 
and the lower edge of the cut as unphysical. 

The goal of our analysis is to investigate the threshold 
behavior of the spectral function in the neighborhood of 
the branch point in a small interval around some tem- 
perature T*, to be specified below. We first restrict the 
discussion to the case |p| = 0, which makes our analysis 
easier to follow. There are now two small quantities char- 
acterizing this situation: e = (T — T*)/T* and 8. [It is 
important to note, that it is M(T) and not M(T*), which 
appears in the definition of 8.] Our basic assumption is 
that the inverse propagator G~ (po) can be expanded 
in terms of these two independent variables. Then to 
quadratic order one obtains 

- G^ipo) « o ie + a 2 e 2 + b x 8 + b 2 S 2 + cSe 

~ : , 1 77 [a l e + a 2 e 2 + bi5 + b 2 S 2 ],{\) 
1 + ce/bi 

where a 2 — a 2 + a\cjb\ was defined, after factoring out 
a temperature dependent finite constant. 

We investigate the behavior of the spectral function, 
which can be derived from JQ| by studying the pole struc- 
ture starting at high temperature. In the concrete case of 
the cr-model the tendency to restore the chiral symmetry 
with increasing temperature reduces the cr-mass below 
2m 7r {T). For such temperature a is necessarily stable if 
the pions are assumed to be stable. 
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When decreasing the temperature the location of the 
CT-pole moves along the physical real axis to the two- 
body decay threshold from below. Let us assume that 
this pole reaches the branch point at T = T* . This is 
the temperature T* around which the behavior of the 
spectral function will be investigated. 

The stability of the O(N) multiplet enforced by the 
chiral symmetry is now postulated, without reference to 
any particular symmetry, for the general discussion. We 
assume that for e > in addition to the approximate 
stability of the particles a\ and a 2 , a is stable. This 
requires a single pole, <5+, for G a on the physical real 
axis below the threshold (1 > 6 > 0), corresponding to a 
stable state. The requirement, that no other positive root 
should be present is fulfilled only if a\ < 0, bi > 0, 62 > 0. 
This choice implies the existence of a root S < on the 
unphysical real axis for both signs of e. This pole stays 
away from the threshold even when e — > 0, see Fig. ^ 

When the temperature is decreased to the region e < 0, 
G^ 1 has no roots on the physical real axis. The pole, 
which did correspond to the stable particle moves con- 
tinuously over the threshold to the unphysical real axis. 
In this temperature range this solution is denoted by 5— . 
This means that for e < one has two poles (S-,5) on 
the unphysical real axis for sufficiently small |e|. 




FIG. 1: Schematic presentation of the pole evolution. 

In this approximation the spectral function p(e, 5) = 
—lmG a (po)/n has the following parameterization above 
the threshold p 2 , > M 2 : 



M 2 
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JT 2 



1)(5- 
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(the expression depends on e also indirectly, through the 
5 ■ 5± term). The pole terms associated with S + and 5- 
in the respective temperature ranges e > and e < can 
be uniquely separated in G a . The approximate spectral 
function based exclusively on this singular "pole" contri- 
bution, above the threshold M reads: 



Psins ~ b 2 (6 ± -6)n X \S\ 2 +5* ± 



(3) 



The maximal threshold enhancement clearly occurs for 
e = 0, when p s - ms ~ 1/| S\, which has been found also 
in Ref. 0- Since 5 + (e) and — 5_(e) in the respective 
temperature ranges are equal to linear order in e, the 
mirror variation to 0(e) of p s i ng around T* is deformed 
only by the temperature dependent first factor on the 
right hand side of Eq. © . 



The threshold enhancement sets in in a universal way, 
reflecting the motion of <5+(<5_) along the physical (un- 
physical) real axis towards the position of the threshold. 
Concerning universality based on the expansion given in 
Eq. Q three further remarks are in order: i) the linear 
dependence of S± on e is universal, which determines the 
leading term of the enhancement according to the second 
factor in Eq. (|3J), ii) when looking at the first correction 
to this leading behavior the one-pole description is not 
sufficient, 5 plays an important role (if the expansion is 
pushed to higher orders, further roots would appear), 
iii) it is obvious that the expansion given in Eq. Q is in- 
dependent of the low temperature quasiparticle spectra. 

Another characteristic temperature value may be de- 
fined, corresponding to the temperature where 5— and 
5 become degenerate on the unphysical real axis, when 
the temperature is gradually decreased further. It is de- 
termined by the smaller negative root of the equation 
b\ — 462(0160 + a 2£o) = 0. For smaller temperature values 
(e < eo) the two poles become complex. One of the poles 
moves into the second Ricmann sheet and can be inter- 
preted as a resonance of finite lifetime if the temperature 
becomes sufficiently low. The other pole is its "mirror" 
located on the continuation of the second Riemann sheet 
through the unphysical real axis into the positive half- 
plane, and does not have any physical meaning. This 
picture fully agrees with the qualitative pole trajectory 
found in p. It could happen that the validity of the 
quadratic expansion of G^ 1 breaks down already before 
the degeneracy of the two poles on the unphysical real 
axis occurs and <5_ may stay on this axis until the tem- 
perature is decreased to zero. This alternative seems to 
be realized in the calculation presented in 0- It would 
be very interesting to see the presence of such a pole in 
the low temperature representation of the pion-pion scat- 
tering amplitude based on chiral perturbation theory . 

The extension of the investigation to the case of 
nonzero spatial momentum |p|/M, regarded as a new 
small parameter can be readily done, by assuming the an- 
alyticity of G„ also on this quantity. Now the threshold 
parameter S takes its original Lorentz invariant form: 5 = 
•yl — p 2 /M 2 , where p 2 is the squared four-momentum. 
At finite temperature, in the propagator also a separate 
dependence on the spatial momentum appears. There- 
fore the expansion Eq. (JJJ will contain also terms of the 
form (p 2 /M 2 ) l e k 5 m (note, that e remains as it was in- 
troduced for |p| = 0!). To lowest order (single pole ap- 
proximation) one keeps terms linear in e and 5 and the 
coefficient d of a new term d x |p| 2 /M 2 has to be com- 
puted. 

The coefficient d determines the shift in the value of 
the temperature for which the maximal threshold en- 
hancement of psing occurs at small finite |p|: T*(|p|) = 

■^t]. Since the term proportional to p 2 is real 



T *[ 1 + 4^7T^ 

even after the continuation above the threshold, |p| and 
e play similar roles in the spectral function. For fixed 
T, it is d, which controls the location and height of the 



3 



maximum of p s i n g when |p| is increasing. For instance, 
one finds S max » l^e + d\p\ 2 /M 2 \/b 1 . 

Below, we illustrate the validity of the simplified repre- 
sentation discussed above for the scalar-isoscalar spectral 
function of the linear a model, obtained to leading order 
in an expansion with respect to the inverse of the num- 
ber of the Goldstone bosons (large N approximation). 
For this approximate solution one checks explicitly the 
existence of the power series of G^ 1 in 6, |p| and e. Fur- 
thermore, the ranges of validity in |e| of the single and 
the two pole descriptions of the spectral function can be 
analyzed by comparing Eq. (J2J and Eq. @ with the com- 
plete expression of the spectral function in the model. 

The first goal in the second part is to sketch the way 
the a propagator given by Eq. (19) of Ref. 6f can be 
cast for po < M = 2m 7r (T) and |p| = in the form 
presented in Eq. (1) of the present paper. Employing an 
obvious notation, the renormalized couplings {m 2 Rl A/?, h) 
of this model were fixed at T = at values which ensure 
m w (0) = 140 MeV, f n = VN<f>(T = 0) = 93 MeV, m a = 
3.95/tt, m a /T a ~ 1.4. The temperature dependence of 
the vacuum expectation value of the O(N) field, $(T), 
was determined by the equation of state (its renormalized 
form appears in Eq. (11) of Ref. |6j): 



$ 2 (T) 



n(u k ) + - 



1 



$(T) 



-0, (4) 



where n is the Bose-Einstein factor and w 2 = k 2 + m 2 . 
The expansion of m^iT) to quadratic order in e is pro- 
vided through the combination of the relation m 2 (T) = 
h/$(T) with the expansion of <f>(T) around <E>(T*). 
$(T*) is found from Eq. (31) of Ref. 0. 

The selfenergy of a is determined to leading large N 
order by a sum of the chain of pion bubbles, characterized 
by the function b > (po) as shown by Eq. (5) of Ref. @ 
(the index > refers to the first Riemann sheet). The 
expression of the bubble above the threshold given in 
Ref. is written, after continuing to the physical real 
axis below the threshold and some convenient integral 
transformations, in a form where the expansion in powers 
of 5 and e actually starts from: 



4TT Z b>(p a ) = - In 



1 m 2 (T) <5arccos(<5) 



M 2 



DC 

/ 



dx 



+ 5 I <lv 

o 



[(l + x 2 )Vl + ( 5 2 x 2 ]- 1 



(5) 



where m 7r (T) and £ = m w (T)/T carry the explicit T de- 
pendence (Mo is the normalization scale of the theory 
defined at T — 0). One performs first the expansion in 
powers of 5 to quadratic order. For the second integral 
in Eq. (J5J one finds to linear order (| — S)n(m 7r (T)). 
Next, the remaining explicit pion mass dependence (also 
in £!) is expanded around T* with the procedure sketched 
above. With these hints one easily establishes the ap- 
proximate expression Eq. Q), obtaining numerically the 



following values for the coefficients: ai = —18.32, 0,2 
-23.14, 61 = 8.59, b 2 = 14.48, c = 9.67. 
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FIG. 2: The spectral function p(po,e)/^ in various approxi- 
mations as a function of po at different reduced temperature 
e, and |p| — 0. 

In Fig. [21 the spectral function and its universal ap- 
proximations are shown for three different temperatures. 
The curves labeled "exact" are calculated from the com- 
plete expression derived in Ref. 6j. They are compared 
with the two approximate forms described above. All 
temperatures were chosen below T* (above it the situa- 
tion is qualitatively similar). For the highest value of e 
the pole S- lies almost on the threshold, for its smallest 
value (e = eo ~ —0.077) one has 5- = 5. When calcu- 
lated from the exact expression one obtains: eo = —0.13. 
According to these curves the region where the one-pole 
description gives a good approximation is rather narrow 
in e and in particular in S. The correction due to the 
second pole contribution is quite important for increas- 
ing frequency values. For e = eg the line corresponding to 
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Psing was omitted, since its expression becomes meaning- 
less. It is, however, remarkable how good approximation 
can be constructed using Eq. 0J. The range of validity 
of the quadratic expansion of G~ x in e includes also the 
degeneracy temperature of the poles, persisting even for 
somewhat lower temperatures when both poles are com- 
plex and lie equidistantly from the threshold. The spe- 
cific form of the spectral function, p ~ |(5|/(|(5| 2 + consi) 2 , 
characteristic at the degeneracy temperature is expected 
to occur even in the exact solution of the TV = 4 linear a 
model with physically relevant couplings. 

For the case of the a— it system, the coefficient d, which 
determines the influence of the explicit |p|-dependence, 
comes uniquely from the temperature dependent part of 
the pion bubble, b^, = b^,{po, |p|). Since the T = part 
of the expression of the bubble is Lorentz invariant, it 
remains unchanged when expressed with the complete 
expression of 6. The expression of b^, for p ^ 0, see e.g. 
Refs. pj |9( can be cast in the form 



b> 



dt n(tm„) 



In 



^ + |p|VF^T 



2,2 



|p|VF=T 



(6) 

For the threshold behavior of the spectral function to the 
lowest (linear) order in (e,6, p 2 ) , one can use 6 = in 
the integral, that is for the quantities scaled by m ff (T) 
one can put: p 2 = 4, p\ = 4 + p 2 . With a suitable trans- 
formation of the variable of integration we can ensure the 
smoothness of the integrand, obtaining: 



b> = 



dt , 1 
In — 



exp [ - f V 1 



4tt 2 £ 



1-expf-A/l 



•IpI 



IpI 



2 coth 2 



ty 



Performing a partial integration one obtains the first two 
terms of the asymptotic expansion of the integral valid 
for small |p|, relevant for our analysis. Next to the un- 
changed |p|-independent term with the complete 5 one 
can calculate the coefficient of p 2 by performing the t- 
integral with the value of £ determined at T*. For the 
actual couplings d = 1.65 is found. 

The assumed expansion of the inverse propagator 



around the threshold in powers of e and S resembles the 
starting point of the Landau theory of second order phase 
transitions. Like in the Landau theory, our purpose has 
not been to prove the analytic property of the relevant 
quantity, but to explore the consequences of its postu- 
lation. We conjecture the universality of the threshold 
behavior in any channel where a particle pole describing 
at high temperature a stable particle might pass through 
the threshold of its lowest two-body decay. Our discus- 
sion implies that for T < T* the pole dominating ana- 
lytically the spectral function cannot correspond to any 
particle. Whether this remains valid down to T = is 
model dependent. 

A possible generalization of Eq. J2J, in which only the 
quadratic dependence of the inverse propagator on 8 is 
assumed, can be easily written down: 



p(e,P,P 2 ) 



Ri(e,p) 



1 



[fa-l + R 2 (e,p)]*+R s (e,p)[^-l} 



( 8 ) 

The functions Ri are real, at T — they assume constant 
values, e is measured from the value T* determined by 
the condition i?2(e = 0, 0) = 0. A sequence of reduced 
temperatures e(p) = (T*(p) - t*)/T* can be defined as 
the solution of the equation i?2(e(p),p) = 0. Obviously 
T*(0) = T* . Asymptotically one expects the "scaling" 
behavior i?2 ~ [ e ~ e (p)] Q > a > 0. The treatment based on 
Eq. leads to a — 1, but in principle a nontrivial non- 
integer scaling exponent can be envisaged. Note, that 
when there is no explicit symmetry breaking the branch 
point tends to the origin if |p| goes to zero and T* = T c 
independently of |p| [lOj. 

The next step will be to discuss the effect of the finite 
pion lifetime. It shifts the branch-point into the lower 
halfplane. Around it there is still an enhancement (sim- 
ilar to the undamped case) but in the spectral function 
varying along the real axis the effect is smeared out par- 
tially, see [lj. Either a complete two-loop (improved) 
perturbative analysis or the evaluation of the next-to- 
leading large N correction to G„ is needed to reliably 
assess the importance of the finite pion width. 
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